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Abstract 

The twofold degenerate Anderson impurity model [1-4] is investigated and the strong 
electronic correlations of d-electrons of impurity ion are taken into account by elaborating 
suitable diagram technique. 

We discuss the properties of the Slater-Kanamori model [2-4] of d-impurity electrons. 
After finding the eigenfunctions and eigenvalues of all 16 local states, we determine the 
local one-particle propagator. Then we construct the perturbation theory around the 
atomic limit of the impurity ion and obtain the Dyson type equation for the renormalized 
one-particle propagator. Diagrammatic theory has been developed and correlation func- 
tion determined. Special diagrammatic approximation was discussed and summation of 
diagram has been considered. 

PACS numbers: 71.27.+a, 71.10.Fd 

1 Introduction 

The theory of strongly correlated electron systems plays a central role in contemporary con- 
densed matter physics. The essence of the problem is the competition between the localization 
tendency originated by the Coulomb repulsion of d electrons and itinerancy tendency arising 
clS db result of hybridization of electron orbitals. 

The orbital degeneracy can be completely eliminated in solid substances but in many of 
them, for example, new superconductors based on Fe and AnCeo materials orbital degeneracy 
is not completely eliminated and orbital effects are important. For instance, orbital degeneracy 
plays essential role in the Mott metal-insulator transition. Here the effects of Hund's rule 
coupling in our orbitally degenerated model are studied with diagrammatic approach. 

We study the influence of the intra-atomic Coulomb interactions of the two electrons with 
opposite spins situated on the same or different orbitals and intra-atomic exchange is analyzed. 

Our investigation is based on the diagram theory elaborated for strongly correlated electron 
systems as in non-degenerated [5-9, 11- 14] and as in twofold degenerated ones [10]. 

The paper has the following structure. In Sec. 1 we describe the twofold degenerate An- 
derson impurity model. The local properties of our model are considered in Sec. 2. The 
perturbation theory around the atomic limit of impurity ion is formulated in Sec. 3. In this 
section we discuss the process of derealization and renormalization of the dynamical quantities. 
In Sec. 4 the simplest irreducible Green's function is calculated. Sec. 5 is devoted to discussion 
of the Mott-Hubbard phase transition. Sec. 6 is devoted to the conclusions. 
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The Anderson Impurity model with twofold orbital degeneracy has the Hamiltonian com- 
posed one part of conduction electrons - one part of interacting localized and strongly correlated 
electrons and of hybridization term between these two parts [1-4]: 



H = H° + H int , (1) 
H° = H C +H^, (2) 
Hi = E^CtC^, (3) 

kla 

H d = Y, edd t d i<r + UY j n n n li + U'n 1 n 2 + I H Y J dtJta d ^'d 2 a' (4) 

l,a I aa' 

+ r H (d+di i d 2i d 2t + h.c), 

H *nt = ^EiV^dtC^ + ViftJ^), (5) 



kla 



where the local Hamiltonian H% is standard Slater-Kanamori [2-4] form, Cg Zcr is conduction 
electron annihilation operator with momentum k, orbital number 1 — 1,2 and spin a = ±l(f, 4)) 
di a operator for localized d electron. Conduction electron oil — th orbital state hybridizes only 
with the local electron of the same orbital state, n\ a = df a di a , n\ — ^2 n^, is matrix element 

a 

of hybridization. U is Coulomb repulsion between the d-electrons in the same orbital state and 
U' - between electrons in different orbital states. Ih is coefficient of the Hund's rule coupling 
and pair hopping terms, ej(fc) is the band dispersion and e d - is impurity ion energy evaluated 
from the chemical potential fi. N is a number of lattice sites. 

In the following we assume that the symmetry of the system is such that exist the relation: 



U' = U -2I H ,I' H = I H . (6) 

The Coulomb interactions are far too large to be treated as perturbation and they must be 
included in H° - zero order Hamiltonian. The hybridization term (5) is considered as the per- 
turbation of the system. In the following the main ideas of the perturbation theory elaborated 
for non-degenerate strongly correlated systems are extended for degenerated systems. Such gen- 
eralization has been discussed, for example, in the case of twofold degenerate Hubbard model 
[10]. As is known the new elements of this perturbation theory of strongly correlated systems 
are the irreducible correlation functions which contain all charge, spin and pairing quantum 
fluctuations. 



2 Local properties 

In the main approximation of the Anderson model one has free conduction and strongly inter- 
acting localized electrons described by the Hamiltonian H . The localized part of the Hamil- 
tonian, Hq, can be diagonalized by using Hubbard transfer operators x mn — \ m )( n \ where \m) 
is eigenvector of operator H% [5]. 

Because orbital quantum number takes two values 1 — 1,2 the total number of local quantum 
states is equal to 16. 

There are the following eigenvectors of operator if J'. The first quantum state |1) is the 
vacuum state |0) with energy E 1 = 0. There are 4 one particle states with spin S = \ and 

S z = 

|2) = dj.|0), |3) = 4rl°>> I 4 ) = d n\°) and I 5 ) = d ti\°)- The energies of all these states are 
E 2 — E 3 — E4 — E 5 — e d . 
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Then there are six states with two particles. Three of them are singlet states with spin 
S — and others 3 triplet states with 5 = 1 and S z = —1, 0, 1, 

| 6 ) = ^(d+ t dj. - 4 t 4J|o), |7> = + 4 t 4)|o), 

|8> = 73«4 - <4 t )|0), |9> = d+d+\0), 

|io) = ^K d n +dJi4)|o>, |ii) = <4Jo). 

The eigenvalues of these quantum states are 

E 6 = 2e d + U - I' H , E 7 = 2e d + U + I' H , E 8 = 2t d + U' + I H , 

Eg = E\q = En = 2e^ + U' — In- 
Then there are four states composed from three particles 
|12) = |0), 1 13> = d^djdjj.|0), 

1 14) = rf+«|0), |15) = 4 t 4<|0), 

with energy value E i2 = Ei 3 = E 14 = E 15 = 3e d + U + 2U' — Ih- 
The last local state is singlet 

1 16} = dJdJ^dJJO) with energy value E w = Ae d + 2U + AU' - 2I H . 
When equalities (6) take place we obtain more simple forms: 

E 6 = E 8 = 2e d + U - I H , E 7 = 2e d + U + I H , E 9 = 2e d + U- 3I H , E 12 = 3e d + 3U - 5I H , 
E 16 = At d + 6U- 10I H . 

The triplet states |9), 1 10} and |11) are the lowest by energy. 

Quantum states enumerated above permit us to organize Hubbard transfer operators x mn 
and establish the relation with fermion impurity operators [10]: 

dt = X 2+l -^ + ^[("l) m X 6 ' 2+ ' +CT + X 7 ' 2+l+ °] + ^X 8 ' 5 - l+ ° + (-!)<+ Y°> 5 -^] (7) 

J^10-<T,5-J-<T _|_ ^_^y+l ax 12+l+er,10+(T _|_ 0^16,15-l+a 

Equation (7) allows to calculate all the local dynamical quantities. For example quantum 
electron number has the form: 

n t = x ^-°W-° + i[ x 6 ' 6 + (-l) l+1 X ^ + (-l) i+1 x 7 ' 6 + X 7 ' 7 ] (8) 

+ 1 [x S,S + + ff( _ 1) I + l x 10,8 + x 10,10] 

10-(T,10-cr _^ ^12+1-0,12+1-0 _|_ ^,12+/+o-,12+/+cr 
_|_ 15— J— <t,15-Z— <r + ^ 16 - 16 

and 

njt _ „ u = _ x 3 + ,,3 + , + (_l)' + l [x 8,10 + x 10,8] + x 9,9 _ ^1,11 + ^U-,,14-1 _ ^6-1,18-1. (Q) 

For r dependent quantity -A(t) = e THo Ae~ rH ° we have the equation: 

n«(r) - nii {r) = X W+l - X W+l + (-1)'+ V' 10 ^"^ (10) 

+ ^10,8 e T(£ 10 -£8)j + ^9,9 _ ^11,11 + X 14-Z,14-J _ ^16-i,16-/_ 

The correlation between quantities with different orbital numbers is determined by the equation 
(1 = 1,2): 

KM - n^r))(n n (0) - n n (0)) = 5 w [x W+l + X W+l + x 14 "'' 14 "' + (11) 

^16-1,16-lj + (_iy+l'^fi e T(Ea-E 10 ) + x 10,10 e r(Bio-E 8 )] + ^9,9 + ^11,11^ 
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EfatOO - n 4 (r))(n n (0) - n n (0)) = -(e"^ 2 + e"^ 2 + 2e"^). (12) 

In special case / = 1, V = 2 we have: 

(m t (r) - n u (r))(n 2t (0) - n 2 +(0)) = -[ X ^ e ^s-E 10 ) + ^0,10^-*)] + ^9,9 + ^ii.n (13) 

which is the d electron susceptibility [2-4] . 

We now define the Matsubara one-particle Green's function of localized d-electrons: 

g°(lar, I'a'r') = g^r - r') = -(T^(r)^(r')) , (14) 

where d la (r) = e TH °d la e- TH ° , d~ la (r) = e THa d^e~ rH °. 

The Fourier components of this Green's function are: 

9 (0) W = k^ (0) K). (15) 
Using (8) and the properties of Hubbard operators we obtain the equation for local function: 

glvo'i^n) = T 771 "FT + - T 771 p~ + ( 16 ) 

1 e -/3£ 2 + e-' 9 ^ 7 3 e - ^ 2 + e _/3B9 e"/ 3 ^ + e -PE 12 
2iu n + E 2 - E 7 + 2 iu n + E 2 - E 9 + iu n + E 6 - E 12 + 

1 e-/ 3 ^ + e -/3£i2 3 e -/3S 9 + e -/3Bi2 e -/9Bi2 + e ~/3S 16 

2 iw n + E 7 - E 12 + 2iu n + E 9 - E 12 + iu n + E 12 - £ 16 ^' 
where Z Q is partition function in atomic limit 

Z = e - pEl + 4e- /3E2 + 2e-^ Ee + e"^ 7 + Se"^ 9 + 4 e - /3£l2 + e"^ 16 . (17) 

The spectral function of impurity d-electron in local approximation is equal to 

A^\E) = -2Img°{E + iS), (18) 

where g°(E + iS) with 5 = +0 is analytical continuation of the Matsubara to retarded Green's 
function. 

Using (14) we obtain 

A^(E) = ^{(e-^ 1 + e- pEa )6(E + E l - E 2 ) + (e'^ 2 + e~^)5(E + E 2 - E 6 ) + 

( e -^6 + e -m^s(E + E 6 - E 12 ) + he~^ + e~^)5{E + E 7 - E 12 ) + 
\{e~^ + e-^)5(E + E 9 - E 12 ) + (e~^ + e~^)5{E + £ 12 - £ 16 ), 

with property 

/ A(°)(£)d£ = 27r. (20) 
J 00 
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3 Derealization processes 



We use the perturbation theory elaborated previously for strongly correlated electron systems 
both of non degenerate [5-9,11-14] and of degenerate forms [10]. We study the process of 
renormalization of Green's function resulting from intra- and inter-orbital flips of tunneling 
electrons. 

The full Matsubara Green's function in the interaction representation for conduction and 
impurity electrons are: 

G(klar\k'l'a'r') = -{TC ka {r)C %IVal {r')Um c o, (21) 
g(lar\l'a'r') = -(Td ltT {r)d lw {r')Umi- 

The anomalous functions are defined as 



F(klar\k'lWr') = -(TC^C^Jr'MP))*, (22) 

F{klar\k'l'a'r') = -(TC^C^Jt'MP))*, 

f(lar\l'a'r') = -(Td la (r)d l/a/ (r')U(p)) c , 

f(lar\lWr') = -{Td la {r)d Va ,{r')Umi- 

Here r and r' stand for imaginary time with < r </?,/? is inverse temperature, T is 
chronological ordering operator. 
The evolution operator is 

f/(/3)=Texp(- f H int (r)dr). (23) 

The statistical averaging is carried out in (21) and (22) with respect to the zero-order density 
matrix of the conduction and impurity electrons. Index c means connected diagrams. 
In the zero order approximation we have 

16 16 

H L = Y. E nX n \ Er = l, (24) 

n=l n=l 

G\% a ,(kk>\r ~ t') = 8 rkl S w S^Gt\k\T-r'), 

r>(P)(u\- ^ 1 (2n + 1)tt 

iuj n - e(k) P 

and g^{iu n ) is determined by the equation (16). 

Hybridization between the conduction and d impurity electrons results in renormalization 
of their propagators. Because the number of conduction electrons N is much larger than the 
single impurity state, the effect of the latter on the conduction band scales as jj. 

The renormalized conduction electron propagator is 

G laVa ,(kk'\\oj n ) = 5^,8 w 5 ac ,GfJ (k\\uj n ) + (25) 
V*V-.., 

-GfJ (k\iu n )gia,l'a' (^n) G£ 0'\^n), 



N 

where gi a i'a>(}u n ) is the full impurity electron propagator. 

A similar equation holds for the anomalous function of conduction electrons in supercon- 
ducting state: 

V*V- 

Fi cVc ,(k,-k'\iu n ) = H N k ' 1 ' GfJ (k\iu n ) fiai'a' i}0J n ) Gfl (-k'\ -iu n ). 



<: 

lar 



:=< 

'crV 



A 



A 



Vi G (o) 



V? 5 



Y 



Vy Gf(0) V^ 2 
>^ a— >: 

1 2 



>: 

lar 



y 



A 



V* G ( ) V 2 f 



+ 



Y 



V x G (o) V 2 * g 



Figure 1: Dyson type equation for Green's function of impurity electrons. A, Y, Y are corre- 
lation functions. 



The equations for the full functions g and / of impurity electrons have the diagrammatical 
form shown in Fig.l. 

The structure representative of the diagrams in Fig. 1 is given by the following equation 



fel k 2 



^rY,\ V k 1 lf G hli(^\ iUJ n)Sl 1 l 2 5a 1 a 2 = Sl^S^G^W 



(0) 



A 



(26) 



where 



A Y io; n - e (k) 



The renormalization quantity is 



1 



Glal'a'i^n) = J2 V kl V i'l' Gl ° l '<?' ( kk ' M ' 
kk' 



(27) 



(28) 



In the Fig. 1 the double dashed lines with arrows depict renormalized g and / propagators 
of localized electrons and solid thin lines represent Q° function of conduction electrons. The 
function V\ means V^ x and summation by repeated indices is assumed. 

A and Y are correlation functions. They contain a sum of strongly connected irreducible 
diagrams. The simplest examples of such diagrams are shown on Fig. 2. 

The analytical form of equations in Fig.l is the following: 
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G^ilfJT, 1\(J\T\\1 2 g 2 t 2 , l'a'r']G^ [hh^^lhhaiTi] 



I a t 




Figure 2: The simplest examples of correlation functions A and Y. 



fi a i' a >(iuj n ) = Y lal , a ,(iu n ) + A il(T1 ; (T (-ia; ri )^ ( i ] i (-ia; n )/ h(7li / (J '( iw n) + 

Yi ahai (i^n)^ (iWn)^iiailV ( iuJ n) • 



(29) 



This system of equations is rather general and admit different phases. We shall discuss one of 
the most simple form with singlet superconductivity on the paramagnetic background. 
For this special case we use the new notations (a = — a): 



gial'a'^n) = <W#" ( iuJ n), fieri' a' i}^n) = 5 a ^ fj a (iu n ) , 



Hal', 
(0)/- 



(30) 



g^ 



By using these definitions we obtain: 

gf(iu; n ) = Af(iu; n )+A^„)^^ (31) 
fli^n) = Yg{^n) + m-^n)S l ^°\-iu; n )f^(i^ + Y^ (iw n )C^ l(0) (iu) n )g l Z (iu n ) . 

In the absence of orbital degeneracy this system of equation has the known solution [14] 



A CT (iu; n ) - gf\-iu n )[k CT {iuj n )K a (-iu n ) + Y^(iu) n )Y St7 (iu n )] 



foofar, 



YrJrriiUn) 



fad(m, 



Y a S\u n ) 



da{iu n ) ' '"' "" d a {iu n ) ' 
d a (iu) n ) = (1 - A CT (iw n )^ 0) (io; n ))(l - A^-iu n )gf' \-iu n )) + 
^ 0) (iw n )^ 0) (-io;„)y^(ia; n )y ff(T (ia; n ). 

Solutions of the equation (31) for the normal state of the degenerate system has the form: 

Ai 1 ^) - GT\iu n )[K^{iu n )Kf{iu n ) - A?(iw„)A»(iu>„)] 



(32) 



9a (^n) 



dJiu n 



9 2 a\iUJ n ) = ^j^, d^U n ) = (1 - Q^(iu; n ) X 

A^(io; n ))(l - gW(iuj n )Al 2 (iu n )) - g^\iu n )gT\iu n )Al\iu n )Al\iuj n ). 



(33) 
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The other two functions are obtained by changing the indexes 1 -H- 2. These equations are of 
Dyson type. They determine Green's functions through correlation functions A = + Z, Y 
and Y ones. The last three can only be given in a form of infinite diagram series, since the 
exact solution does not exist. 

An example of efficient summation of diagram and determination of the correlation function 
Z, Y and Y is presented on the Fig. 3. 



ham 




v 2 

l 2 a 2 T 2 

Gf iT \lar, ham\l 2 a 2 r 2 , l'a'T']V^G[k 2 ha 2 T 2 \kik^iri]V2 



I'a'r' 



ham 



YfAr-r') 




V 2 

l 2 a 2 r 2 

G 2 ^ ir [laT, l'a l T J \ham,l 2 a 2 T2]V\F[kiham\ — k 2 l 2 a 2 T 2 ]V 2 



I'a'r 



ham 




l 2 a 2 r 2 

G { 2 0)ir [ham, l 2 a 2 r 2 \lar, I'a'r'jV,* Fi-hhamlhh^M 



Figure 3: The main approximation for the correlation functions. The solid double lines with 
arrows depict the full Green's functions of conduction electrons. The rectangles depict the 
irreducible Green's functions of the impurity electrons. 



The diagrams of Fig. 3 differ from the ones of Fig. 2 by the presence of the full conduction 
electron Green's function instead of the bare one of Fig. 2. This difference is the result of 
ladder summation of main diagrams. 



4 Correlation functions 

The simplest correlation function is determined as 

Gr[l,2|3,4] = ^ } (1,2|3,4) -^ 0) (1|4)^ 0) (2|3) +^ 0) (1|3)^ 0) (2|4), (34) 
^ 0) (1,2|3,4) = (Td 1 d 2 d 3 d i ) , <7i 0) (l|4) = -(Td 1 d A ) ,l = fr^n), 

with two- and one-particle bare Green's functions of localized electrons. 

Because the presence of the Coulomb interactions in zero order Hamiltonian, equation (34) 
is different of zero and contains charge, spin and pairing fluctuations. 

The two-particle Green's function #3 is the sum of 4! terms of different time ordered 
electron operators products. The statistical averages of these quantities are calculated by using 
Hubbard transfer operators representation. 
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We need the Fourier representation of these functions 



G t 2 r [lia 1 T 1 ;l 2 (T2T2\hcr3T 3 ;Ucr4n} = ~L G 2 rr [ha^u^ k^i^lk^ius] l 4 a^u 4 } x 

P U11UJ2LU3L04 

- iui 1 Tl — iW2 T2 +IOJ3 T3 + I0J4 T4 



e 



gf\k^n\l 2 a 2 r 2 ) = ^g?\l 1 ar,l2<J2\iuj 1 )e- i ^-^ (35) 

P UJl 



(0)/l 7 j • X r r /. X ^W2^01CT 2 / ' 

01 (tio-i; «2^2Ni) ~ ^i«2^i ( T2^( 1 wi) = — o — • — — et + • 

Z VlWi + hio — hi a 1 



VjJ\ + i?2 — Eg iu± + i?g — -&L2 

G 2 rr [/icriia;i; ^o^ic^ I ^303^3; = g^ihviiwi, l 2 (T 2 ioj 2 \l 3 a 3 ioj 3 ; l A a A \u A ] - 

f38(ui + lo 2 - u 3 - u)4)[f38(u;i - co^g^ihai, h^iu^g^ {1 2 g 2 \ l 3 a 3 \\uj 2 ) - (36) 

(35(ui - u 3 )gf- ) (ha 1 ; l 3 cr 3 \ioj 1 )gf > \l 2 a 2 ; ha 4 \iu 2 )}. 

There exists the law of frequency conservation 

GY[hoiiwi\l 2 V2iw2\hazVjJz\haAVjJ^ = (37) 
(35(ui + uj 2 — C03 — uj^G^iha^uji, l 2 a 2 iu 2 \l 3 a 3 iu 3 ; Ua^u;^. 

The statistical averages of chronologically ordered products of the electron operators of the 
function g^ have different weights of the form 6 ^ n , where E n are the energies determined in 
previous section. Because Eg is the lowest energy the weight e~^ E9 is the main of them and 
only such terms are taken into account. 

Just such considerations determined us to use instead initial exact equation (16) for zero 
order Green's function g^ V(T i the approximate value (35). Zero order partition function Z (17) 
concomitant is approximated as 3e -/3£;9 . 

For example the contribution to function g^ [ii<7iiu;i; l 2 a 2 iuj 2 \l 3 a 3 \uj 3 \ ^o^iu^] with time order 
j3 > Ti > r 3 > r 2 > r 4 > and with weight e - ^ 9 is 

^li3^2u(^ ( ^0'l(T3^<T2IJ4 ^0\,—OZ^Ol,—Olfi<Jl03 ^<Tl (T3 ^2(74 ^fl^^^fs) ^1324 
(^3-Zi-« 3 ,0^3-i 2 -/4,0 + (-l)' 1+U ^i/ 3 ^2u)(^ cr l (7 4^<T 1( T3^a2 ( T4 + " CT 3 ^2 - (74^1(74) -^1524 _ (38) 

{-^) h+l4 h-h-hflh-h-hfl{^ a l a ^o 1 a: i 8cT2<T4 + ~ CT 3 ^2,-CT 4 ^(71 £74)^1324' 

where 



/9 ri r 3 r2 



J l324 = / rfr i/ rfr s/ dr4e (£9 - Bl2)(n+r2 - r3 - r4) e iWiri+iW2T2 - itU;ir3 - iW4r4 , 



/g_ = 1^ y d Tl Jrf T3 Jrf T2 J rfT4e (B9-fii2)(r 1 -r4)+(E6-Ei2)(r2-r 3 ) e ia; 1 r 1 +i W 2r2-ia;37 5 -ia;4r4 ) (39) 


l[f 2i = <—^-L J dn jdr 3 jdT 2 jrf T4e ( S 9- E l2)(n-r 4 ) + (i?7-Sl2)(r2-r3) e ic. 1 r 1 +i W 2r 2 -ic,3r3-^4r4 > 


These 4-fold multiple integrals by time variable r can be transformed in contour integral 
by using the method of Claude Bloch [15]. With this purpose it is necessary to introduce the 
exponential form 

e (l3-T 1 )Eo + (T 1 -T- i )E 1 +(T-i-r2)E2 + (T2-T4)E- i + {T4-0)E4 
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which must be compared with exponential form of our integrals I^ 24 - Comparison with I^ 24 
give us the result 

Eq = —Eg, E2 
r(i) 



= -E 9 + iui - iu 3 , E 4 = -E 9 + ifi, Ei = -E 12 + icoi, 
E 3 = -E12 + iui + iu 2 - iu 3 , ^ = ^1 + w 2 - u 3 - u 4 . 



(41) 



Our integral IW 24 ls transformed in the contour integral 



jd) = J_ J_ 



2niZ J (z + E )(z + E 1 )(z + E 2 )(z + E 3 )(z + E 4 ) ' 



(42) 



where contour C + surrounds the real axis in the positive direction. The integrals 1^ and 
have the same form (42) but differ in the definition of energy E 2 . For 1^ the energy 
E 2 = — Eg + \lo\ — iu 3 and for E 2 = -E 7 + iu>i — icu 3 . Other parameters coincide. 

The contour integral (42) is evaluated by the method of residues. The simple results are 
obtained when the parameters E n are different. The existence of multiple poles is possible for 
the special values of frequencies u> n . 

For example in the case when oj\ — u 3 = and fi = 0we have E = E 2 = E 4 and the pole 
z = —E is 3-fold multiple with the residue 



2 \(z + E^z + E 3 ) , 



(43) 



z=-E 



To find all possible multiple poles we consider different values of frequencies using the 
identity 1 = 5(u) + ip(oo), where i/j(u) = 1 — S(oo). For example we consider the possibility when 
fl can be equal to zero and u>i — ui 3 . We have the identity: 



1 = (8(Q) + V(n))(<J( Wl - u 3 ) + ^ - u 3 )) = 
<y(fi)<J(wi - co 3 ) + 5(n)^(u 1 - co 3 )) + ^(fi)<J(wi - cos) + ^(^)V'K - w 3 )). 



(44) 



The first term in the right-hand part of this equation admits the existence of triple pole, the 
next two terms admit double poles and last term admit double and single poles. 

We shall take into account these residues, statistical weights of which is 6 z 9 , and shall 
omit the other ones. In such approximation we have 



Zal 



(i) 



1 



-pz 



0^1324 = ^( )^1- W 3) 



8(Q)i/j(ui - u 3 ) 
S(ui - u 3 )ip(Q) 



e -pz 


(z + E 1 )(z + E 2 )(z 


+ E 3 ) 


e -pz 




(z + E 1 )(z + E 3 )(z 


+ E 4 ) 


e -pz 




(z + E )(z + E{){z 


+ E 3 ) 



(z + EJiz + Es) 
P -pz \ 



z=-E 



+ 



-E 



+ 



-Eq 



z=-E 2 

ll){Q)i>(u-y - UJ 3 )lp(u 2 - U 4 )[ 
v-PZ \ 



+ 



(z + E Y(z + E 1 )(z 


+ E 3 ) 


e -pz 




(z + E Y(z + E 1 )(z 


+ E 3 ) 


V'WK - 


u 3 )5( 


e -f>z 




(z + E 2 y(z + E 1 )(z 


+ E 3 ) 


e -pz 





z=-E 2 



z=-E 4 l 



z=-E l 



{z + E )(z + E 1 )(z + E 3 )(z + E 4 ) 



+ 



-E 2 



z + E 1 )(z + E 2 )(z + E 3 )(z + E 4 ) 
e -pz 

Jz + E )(z + EJiz + E 2 )(z + E 3 ) 



+ 

+ 

+ 
x 
+ 

+ 



z=-E 



-E 4 
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The contribution of other poles is negligible. Our next approximation consists in preserving, 
in the case of low temperature, of the main part of the second derivative (43) just of the form 



AI 



(46) 



2Z (E 1 — E )(E 3 — E ) 

This contribution together with contribution (36) of the product of one-particle Green's 
functions determines the main part of the correlation function. This part is designed as Q^ irr . 

After some transformation and summation of different contributions we obtain the main 
approximation for the correlation function: 



G 2 r [licr 1 iu) 1 ; l 2 a 2 iu 2 \l 3 a 3 iu 3 ; ha^iu A } = ^-5(ui + u 2 - w 3 - w 4 )p(iwi)p(iw 2 ) x 

6 

{fi5{ui - oj 4 )S hh Si 2 i 3 [25 



P 



with 



piiuj) 



iu + E 2 — Eq iu + Eq — E 
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(47) 



(4£ 



5 Mot t- Hubbard phase transition 

As has been mentioned above, one example of efficient summation of diagrams which determine 
correlation function Z and A is presented on the Fig. 4. It has the form 



k{x\x') 



-1 




-1 




Figure 4: The main equation for the function A(x\x'). Here x is (/, a, iu). The thin dashed line 
represents the bare local one-particle Green's function and the double dashed the renormalized 



one 



The thin solid line represents the conduction propagator. 



First of all we shall discuss the approximation with zero order correlation function G 2 
Using the result (47) we obtain 

^^GfltoilKTMlI - .i/-/;,..i*(o) (it. i \ X x .x...r^ ; ,..M2„(o) 



x a 1 \uj 1 -,l'a'iuj]g\ll l {iuj l ) = ~5 aa/ 5u>[p(iu)] 2 g^\iu), 



(O)irr 



(49) 



G 2 [laiu;haiiui 



^p{iu)p{iui){f35(u - u))6 u > O^.-^SO^O 

h 

-<W Q^-a i}^l)9h,-ah-ai}^l) Q^- a fal ) ) 

-(35(u - u) l )[25 aa >gj;'z r7 (iu) l )gt _ aV „ a (iui)Gf^ a (iui) + 
-<W Gta (iui)gi*i><T (iwi ) Qfl (iu x ) - <5 ff _ CT ' (iwi) </w^ (iwi ) Gf}, (i^i )] } . 



(50) 
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We keep the terms which preserve the spin and have the form 8 aa/ and omit the terms with 
spin-flipp of the form 5 a /- a and also omit the terms which are reciprocally subtracted and differ 
only by the sign of spin. We take into account that the function Qi°\iui) doesn't depend of 
spin index and 

E^Sm^- ( 51 ) 

CI 

As a result of such simplifications we obtain 

G£ 0)OT [ZOTw;ZiCTiiwi|Z 2 a2^ = ( 52 ) 

~0S(u - Wi)p(ia;)p(ia;i)5 CTCT /^ ) (ia;i)^/ CT (ia;i)^° ) (ia;i), 

= <^<WHM + ^-gi°\iu)] + ^(^(iwWH, (53) 
with the following realizations 

An(iw) = mi(ia;)+ - g| + - [gi >)] fln(ia;), 

A22M = m.M + ^^f (i w ) + ^M[^°)(i w )] 2 fe (ia;), (54) 

A 12 (iw) = ^|^^ 0) (i^)^ 0) (ia;) 5l2 (ia;). 
We take into account the Dyson type equation 

An(iw) -^ 0) (iw)(An(ia;)A 22 (ia;) - A 12 (iu;)A 21 (iu;)) A 12 (iw) 

" df>) " "dM"' (55) 

= (1 - An(i^)^ 0) (^))(l " A 22 (i^)^ 0) (i W )) - ^ 0) (i^)^ 0) (ia;)A 12 (ia;)A 21 (i^)). 

We make some generalization by considering function m(iu) dependent on orbital quantum 
number I even if it is really not. The function g 22 is obtained from equation (55) by changing 
indices 1 and 2. 

We have found two solutions of equations (54) and (55). 

The first of them is 

An(ia;) = rai(iw), gu(iu>) = 



G\ M G\ M V ^r 1 ^ (Mr 

a 12 m = a 21 m = ±— =-, g 12 M = ± (0) (0) . 



with the condition that 



mo),. : - m i( iw ) = r (o) r : ~ m ^ luJ ) ( 57 ) 

The second solution is obtained from (56) by changing indices 1 and 2. 

The analytical continuation of obtained solutions in upper semi-plane gives us the possibility 
to determine spectral function of localized electrons 

p u ,(E) = -2Img u ,(E + i5). (58) 
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For example intraorbital contribution has a form 

Pn(E) = _ 2 "f^ + 'l , (59) 

where 

Img[°\E + id) = -7rpo(£)N 2 . (60) 

The quantity pn(E) differs from zero thanks the existence of the matrix element of hy- 
bridization and of the zero order density of states po(E). For E — po(0) is positive and state 
of the system is metallic. 

Interorbital contribution to the phase transition is determined by the value 



p 12 (E) = -Img 12 (E + i5) = : ^ ^ (61) 



2(E + AE 1 )(AE 2 - E) 
IniGf\<S)ImG { ?\<S){AE 1 + AE 2 ) 

where 



AE\ = E 2 — Eg > 0, AE 2 = E\ 2 — Eg > 0. 
This quantity is positive for —AE 1 < E < AE 2 . 

For these energy values the state of the system is metallic. The appearance of spectral 
weight et the Fermi level is considered as a definition of Mott transition. 



6 Conclusions 

Diagram approach for investigation the properties of twofold degenerate Anderson impurity 
model has been elaborated. 

First of all the eigenfunctions and eigenvalues of energy of the localized d— electrons part of 
the Hamiltonian have been determined. Their dependence of intra and inter orbital Coulomb 
interactions and of Hund rule coupling constant was established. 

Perturbation theory around the atomic limit has been developed and Matsubara Green's 
functions as in normal and in superconducting states has been defined. 

Dyson-type equations for these functions have been established for both states but detailed 
solutions were discussed only for normal state supposing additional investigation in the next 
paper. 

Because the main elements of our diagram technique are the irreducible Green's functions 
we have undertaken the determination of simplest two-particle irreducible Green's function and 
determined its dependence of the spin and orbital quantum numbers. This quantity has been 
determined only in the low temperature limit. 

Having this quantity and summing some class of diagrams we have obtained the h-icVa 1 
correlation function. 

We found two solutions for the renormalized Green's functions of the d— electrons and 
determined the spectral weight. 
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